Oblique solitons in the flow of polariton condensate past an obstacle 
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Formation of oblique solitons by a flow of polariton condensate past an obstacle is considered. The 
flow is non-uniform due to a finite life-time of polaritons what changes drastically the conditions of 
formation of oblique solitons compared with the atomic condensate case. The theory shows that the 
polariton solitons can be generated by a subsonic flow in agreement with the recent experiment [A. 
Amo et at, Science, 332, 1167 (2011)]. Geometric form of oblique solitons and other their parameters 
are calculated analytically and analytical results are confirmed by numerical simulations. 
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1. Introduction. Semiconductor microcavities, in 
which strong coupling of excitons to photon modes leads 
to the formation of microcavity polaritons, have become a 
focus for the study of polariton condensation and related 
physical effects in solid state systems 0, Q • Polaritons 
possess an extremely small effective mass m po i of the or- 
der of 10 -4 that of electron which allows for their conden- 
sation at temperatures about a few kelvins and more. All 
parameters of the condensate can be tuned with the use 
of resonant lasers and this advantage of polariton conden- 
sate compared with the atomic case has been exploited 
in the recent observations of macroscopic coherence |3j- 
01, quantized vortices superfluid flow past an obstacle 
[7H10|, and persistent superfluid currents [llj . 

Theoretically, different regimes of the flow of atomic 
condensate past an obstacle in two-dimensional (2D) ge- 
ometry can be distinguished depending on the Mach 
number M = uq/c s (uq being the velocity of the inci- 
dent uniform condensate and c s the sound velocity in it): 
(i) for < M < 0.43 the flow is superfluid and no excita- 
tions are generated 12, 13j]; (ii) for 0.43 < M < 1 vortices 
are generated by the flow what corresponds to the usual 
mechanism of loss of superfluidity; (iii) for 1 < M < 1.44 
the channel of Cherenkov radiation of sound waves opens 
which leads to formation of the interference wave pattern 
("ship waves") outside the Mach cone vortices 
are still generated and they are located inside the Mach 
cone forming oblique vortex streets; (iv) for M > 1.44 
the oblique vortex streets are transformed into oblique 
dark solitons 17J which are surprisingly stable and this 
is explained as a transition from absolute instability of 
dark solitons to their convective instability flil - l20l |. 

The dissipative nature of polaritons changes essentially 
the properties of wave patterns generated by the flow of 
polariton condensate past an obstacle. For example, the 
amplitude of the Cherenkov wave pattern decays very 
fast with growth of the distance from the obstacle even 



in the pumping regime as it was predicted in [2l| and 
observed experimentally in Q. Situation changes even 
more drastically in the flow of the polariton condensate 
outside the pumping region [§] or in the free flow of the 
condensate cloud [uj. As was found in {H, the stable 



oblique solitons are generates by a subsonic flow with 
the Mach number M ~ 0.6 and the above mentioned 
region (ii) becomes extremely narrow in the polariton 
condensate flow. We note that in this experiment the 
condensate flow was strongly non-uniform — its density 
decreased with the distance from the pumping region be- 
cause of dissipation effects and the flow velocity increased 
due to action of pressure, thus the Mach number M be- 
came a function of the space coordinate. In this Letter we 
develop the theory of oblique dark solitons in such non- 
uniform situations and show that the dissipation effects 
relax the condition for convective instability of oblique 
solitons. Thus, for subsonic incident velocity and strong 
enough damping, the unfavorable for stability region be- 
comes so narrow that it can be spanned by the oblique 
soliton wave and this explains the results of Ref . . 

2. Theoretical model. Polaritons have a finite life-time, 
and to maintain their steady-state population a contin- 
uous pumping is required. Above a threshold pumping 
strength an accumulation of low energy polaritons is ac- 
companied by a significant increase of spatial coherence 
that extends over the entire cloud of polaritons 0, Q 
which, hence, with account of repulsive interaction of 
polaritons caused by their exciton component, can be 
described by the Gross-Pitaevskii (GP) equation for the 
polariton condensate wave function tp. Due to dissipa- 
tion the polaritons disappear from the system, and these 
losses can be modeled by the effective term dtip = —J"4>- 
As a result, outside the pumping region, we arrive at the 
following generalized GP equation 



iipt + h^ip ~ |V>| 2 V> = U(r)ip - ryV>, 



(1) 



(written in standard non-dimensional units) where U(r) 
denotes the potential of the obstacle and A is the 2D 
(r = (x,y)) Laplace operator. In accordance with the 
experiment [9j , we assume that the polariton condensate 
is created in the region x < with the density po and 
it is put in motion with the flow velocity uo- Hence, 
in the region x > it evolves according to Eq. ([l) and 
satisfies the boundary condition ip\ x — = \/po exp(zuox). 
Numerical solution of this problem is illustrated in Fig. 1. 
Similar pattern was found in numerical solution of a more 
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FIG. 1: (color on line) Distribution of the density plot of 
polariton condensate in the region x > for the value of 
the damping coefficient 7 = 0.02. The obstacle is located 
at the point x = 1.0, y = and modeled by the poten- 
tial U corresponding to an impenetrable disk with diame- 
ter equal to unity. The parameters of the incident flow are 
po = 1.0, mo = 0.65. White curves show the locations of 
the solitons according to the analytical formulae (see (|5J) and 

mi)- 



general model [22j. However, simpler model ([1} is enough 
for our aim to understand the behavior of oblique solitons 
with account of dissipation and it permits us to obtain 
more detailed analytical results, so we will study it here. 
As we see in Fig. 1, the distribution of density decays 
with growth of x and the oblique solitons are created by 
a subsonic incident flow. Due to non-uniformity of the 
flow, the oblique solitons are curved. Thus, our task is 
to develop the theory of these effects and to give sim- 
ple estimates of the main parameters of the flow and of 
solitons. 

It is convenient to transform Eq. ([1]) to the hydrody- 
namic form by means of the substitution 



V> = \fp exp i^i J u(r',t)dr' 



which yields 



p t + V • (pu) = -27p, 



(2) 



(3) 



u t + (u • V)u + Vp + V 



Then the boundary conditions reduce to 

p = po? u = uq at x = 0. 



(5) 



Equation (JXJ) with 7 = 0, U(r) = and uniform back- 
ground density p = p b has the dark soliton solution which 



in the reference system with the quiescent condensate 
(uq = 0) can be written in the form 



p s (x,t) = p b < 1 



1 - V 2 / Pb 



cosh 2 y Pb - V 2 (x - Vt)] 



(6) 



u s (x, t) = V(l - Pb /p s (x 1 1)) 



where x axis is the direction of propagation. The oblique 
solitons of Ref. [17| far enough from the obstacle can be 
considered as solutions ^ transformed to the reference 
system with the obstacle at rest and such a value of the 
velocity V that it is canceled by the component of the 
flow velocity normal to the soliton. We have to find how 
this solution deforms in a non-uniform condensate. 

3. Hydraulic approximation. If the dependence of the 
background steady flow on x coordinate is not too fast, 
then it can be found in the hydraulic approximation with 
neglected dispersion terms in Eq. Q and the time deriva- 
tives put equal to zero. The flow becomes smooth enough 
at some x — x\ where p = p\ and u = u\. The system 
can be easily reduced to a single equation 



(7) 



which solution with the initial condition p — p\ at x = x\ 
gives the background density p b {x) for x> x\. It is given 
in an implicit form by the equation 



x = X\ 



(^pl/ 7 )X(p b / Pl ), 



where 



X(p) = §(Mi - [Ml + 2(1 - p)] 1 ' 2 ) + {M 2 + 2) 1 ' 2 



x In 



sTp- 



(M 2 + 2) 1 ' 2 - Mi 



(M 2 + 2)V2 _ [M 2 + 2(1 - p]V2 
Mi = i*i / \fp~i and 



(8) 



(9) 



u b {p b ) = \Ju\ + 2(p 1 - p b ) 



(10) 



is the background flow velocity. These formulae deter- 
mine the distributions of the density p b and the flow ve- 
locity u b as functions of x. If the incident flow is super- 
sonic, u > y/po, then they are applicable in the entire 
region x > 0, so that x± can be put equal to zero and, 
hence, p\ = po and u\ = uq. For a subsonic incident flow 
velocity uq < ^/po the profiles change very fast near the 
input plane x = and the dispersion effects cannot be 
neglected. However, at the distance about one healing 
length the profiles become smooth enough for applicabil- 
ity of the hydraulic approximation (|9| , (fTOl) . hence we can 
choose the reference point as Xi ~ 1. As we see, the pa- 
rameters of the flow change considerably at the distance 
about y/pxj^i that is of the sound velocity multiplied by 
the time-life of polaritons. Accuracy of the formulae ob- 
tained is illustrated in Fig. 2 by the plot of a local Mach 
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velocity, that is 



FIG. 2: (color on line) The local Mach number M(x) = 
Ub(x) I \J pb{x) as a function of x. The parameters of the inci- 
dent flow are po = 1.0, wo = 0.65; 7 = 0.02. The continuous 
line corresponds to the analytical formulae and the dots to 
the numerical solution of Eq. (0. 



number M(x) = Ub(x)/ \J pb(x) as a function of x. As we 
see, for typical values of the parameters the Mach num- 
ber exceed the critical value M = 1.44 of transition from 
absolute instability of oblique solitons to their convective 
instability very close to the obstacle at the distance much 
less than the size I ~ y/po/j of the condensate and the 
hydraulic approximation is applicable practically every- 
where. 

3. Dark soliton on the slowly dissipating background. 
Motion of dark solitons along non- uniform background, 
e.g., along the ID Thomas- Fermi profile formed by the 
stationary trap with a harmonic pot ential, was studied 
in 23 1 and it was shown in [24I |25| that dynamics of 



dark solitons is determined by conservation of the soli- 
ton's energy. In our dissipative case neither the num- 
ber of polaritons, nor the energy is preserved and the 
stationary profiles ([§)), (fTUj) are created by decaying non- 
uniform flow rather than by the external potential. We 
suppose that if the profiles of the density and the flow 
velocity are varying slowly enough, then the dark soli- 
ton can adjust its local parameters in such a way that 
its form © is preserved with varying density pb = Pfe(r), 
velocity V = V\ \J pb{r)/pi, and x — Vt replaced by the 
coordinate £ normal to the curve y = y s (x) of the soli- 
ton location (determined, e.g., as a curve along which 
the local density changes most slowly). Here pi and u\ 
are taken at the reference point ri which can be chosen 
closely enough to the obstacle and then we can assume 
with good enough accuracy that V\ — u\ sin a where 2a 
is the opening angle of two oblique solitons generated by 
the obstacle (see Fig. 1). The curve y = y s (x) is deter- 
mined by the condition that the local soliton velocity V 
is canceled by the normal component of the local flow 



dy s 
dx 



= tan I 



sm 




Pb{x) 



pi[u\ + 2(pi 



Pb(x))} 
(11) 

It is easy to see that the solution can be reduced with the 
use of Eq. ((U to a quadrature which in a practically useful 
limit V\ <C y/pi can be expressed in terms of elementary 
functions 



y s = Vi + (y l /2 1 )Y{p b / Pl ) 



(12) 



where 



Y(p) 



x In 




(13) 



Then the functions x = x\ + (^/pi/j)X(pb/ pi), y = 
±[?/i + (V\/2"/)Y(pb/pi)] determine in parametric form 
the locations of the dark solitons shown in Fig. 1 by 
white lines. Observed agreement of analytical formulae 
with the numerical results justifies our supposition that 
the parameters of oblique soliton adjust their values to 
the local parameters of the flow. 

4- Stability of oblique solitons in a non-uniform back- 
ground. As is known, 2D dark solitons propagating 
along quiescent background are unstable with respect to 
"snake" instability [261428 1 which leads to decay of dark 
solitons into vortices. However, dark solitons in BEC can 
be stabilized by the flow of the condensate which convects 
the unstable perturbations away from the obstacle 18 1. 
As was shown in [l9j . such a transition from absolute 
instability of dark solitons to their convective instabil- 
ity corresponds to the condition that the flow velocity 
along the soliton exceeds the velocity of propagation of 
the front of the instability region and this front veloc- 
ity equals to the minimal group velocity of the waves 
propagating along the soliton. To apply this theory to 
our problem, we have to generalize it to the case of dark 
solitons propagating through a non-uniform background. 

The spectrum uj(p) of harmonic waves oc exp[i(prj— cot)] 
propagating along dark solitons (here r\ is a coordinate 
along the soliton) was calculated in 28j where it was 



shown that for wave numbers < p < p c , where 



1/2 



p c {V)= 2^/V 4 - V 2 + 1 - (1 + V 2 ) 



(14) 



V being velocity of the soliton moving with respect to the 
quiescent condensate with a uniform density pb = 1, the 
spectrum is pure imaginary, i.e. it corresponds to unsta- 
ble modes, and for p > p c it is real, i.e. it corresponds to 
propagating modes. The analytical expressions for u)(p) 
are known in two limits: for p « p c (see (2^|) and for 
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\p — p c \ < Pc (see 29]). We will use here an accurate 
enough interpolation formula 



3pciv) ^+^^Jy), (15) 



p(V)=b(l-V 2 )/ Pc (V), 

where b = 0.2369 is a fitting parameter. Then the mini- 
mal group velocity equal to the unstable front velocity is 
given by the formula 

v g (V)= 30b (J 1 + f k(1 + k) - (1 + 2k) J 
x ^(1 + + f k(1 + k)- 1 + |k(1 + rc)J 



(16) 



where k = k(V) = 6(1 - V 2 ). 

In case of slightly non-uniform background we get the 
following approximate expression for the velocity of the 
instability front as a function of x coordinate of the soli- 
ton location point: 



v g (x) = \Jp b {x) VgiVx/y/p^), 



(17) 



where the parameters V\ — uisina and p\ refer to the 
reference point near the obstacle at which the hydraulic 
approximation is applicable with desirable accuracy. If 
the component of the background flow velocity along the 
soliton exceeds v g (x), then at this point x the instability 
front is convected away from the obstacle and the oblique 
soliton increases its length. This condition yields the 
inequality 



Pb(x) < pi 



Ml 



V 1 2 /p 1 + 2 + v g -(V 1 /^pl) 



(18) 



which determines the region of x where the oblique soli- 
ton becomes only convectively unstable, i.e. effectively 
stable. The right-hand side of this inequality has a 
minimum for the most deep solitons with V = and, 
hence, v g (0) = 1.44. This gives the upper estimate 
for the x cr such that (fT5| is satisfied for x > x cr : 
Pb{x cr ) — 0-2A5pi(M 2 + 2), and substitution of this value 
into Eq. ([9J yields the expression for x cr , 




F(M0, 



(19) 



1.1 1.2 1.3 1.4 

FIG. 3: Plot of the function F(M) (see Eq. (50 



remarkable, this function takes very small maximal value 
about F max = 0.055. Just due to this smallness the crit- 
ical value x cr is much less than the characteristic size 
of the decaying condensate cloud I ~ y^Po/l- For sub- 
sonic incident flow the distance x cr can be calculated nu- 
merically only but the result is approximately the same. 
Thus, we arrive at the following scenario of formation 
of oblique solitons by the flow of polariton condensate. 
There are two characteristic distances — the condensate 
size I and x cr ~ 0.055 • I differing by more than order 
of magnitude. If x cr is about a few healing lengths at 
the location of the obstacle, then the shadow "seeds" of 
oblique solitons can easily reach the region of the effec- 
tive stability x > x cr and be observed in a wide region 

X cp X I • 

Let us make estimates for the parameters of the ex- 
periment There the sound velocity at the location 
of the obstacle was estimated as c s = 3.5/im/ps and the 
polariton life-time as 15ps which gives the condensate's 
size I = 50/zm as was observed in the experiment. Then 
we get x cr — 0.055 • / = 3/im which is about the observed 
width of the oblique soliton or the healing length. Thus 
the oblique solitons are formed actually in the region of 
their effective stability what explains their formation at 
subsonic Mach numbers of the incident flow. 
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A. Bramati, I. Carusotto, C. Ciuti, E. Giacobino, 
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superfluidity in the cavity polariton physics. This work 
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where 



F(M) = ^ (m - v/0.509(M 2 + 2) 
+ \JM 2 + 2 In [l.729(v / M 2 + 2 - M) 



(20) 



The plot of the function F(M) is shown in Fig. 3. It 
vanishes at M = 1.44, as it should be, and, what is most 
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